Subharmonics and Aperiodicity in Hysteresis Loops 
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We show that it is possible to have hysteretic behavior for magnets that does not form simple 
closed loops in steady state, but must cycle multiple times before returning to its initial state. We 
show this by studying the zero-temperature dynamics of the 3d Edwards Anderson spin glass. The 
specific multiple varies from system to system and is often quite large and increases with system 
size. The last result suggests that the magnetization could be aperiodic in the large system limit for 
some realizations of randomness. It should be possible to observe this phenomena in low-temperature 
experiments. 
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Hysteresis in magnetic systems Q is the prototype for 
understandinghistory dependent behavior in all complex 
systems, pi ESQ] Although the usual example of this is for 
ferromagnets, it can exist in other magnetic systems such 
as anti-ferromagnets^ or spin glasses. 00 

For magnetic systems, when one cycles adiabatically 
between a large positive and negative field H, the mag- 
netization M forms a closed loop in the H — M plane 
because of saturation. There is an implicit expectation 
that even if the field were cycled repeatedly between two 
moderate values, the magnetization would settle down to 
a steady state closed loop, after some initial transient re- 
sponse. In this letter, we show that this expectation is 
often wrong. Many cycles of the external field are neces- 
sary to bring the magnetization back to its initial value, 
and that as the system size grows, so does the number of 
cycles needed. 

We do this for Ising spin glass models: both short range 
and long range models violate the usual assumption of 
simple closed loops. Even in steady state, the magneti- 
zation returns to its initial value only after multiple cycles 
of the field. Another way of looking at this phenomenon, 
probably more relevant to real experiments, is that with 
a steady state oscillatory field the magnetization shows a 
subharmonic component. The number of cycles it takes 
for the magnetization to return varies from system to 
system, depending on the realization of the randomness, 
and can be very large. Our estimates suggest that this 
effect should be observable in nano-particles of spin glass 
material at very low temperatures. 

This result might appear to be rather surprising as 
there are many cases when the strong result of "return 
point memory" 0, [|J can be shown: if a system starts at 
a field Hi and ends at a field H2, its final state is in- 
dependent of the time dependence of H as long as H is 
never outside the interval [Hi, H%\. In the case above this 
would imply at most one transient cycle before the steady 
state closed loop was reached. An example of such a sys- 
tem is a random field Ising ferromagnet where an elegant 
proof of return point memory was given. 8] However the 
proof does not apply to cases where antiferromagnetic 
interactions are also present, such as a spin glass. We 
show below that in this case there is a much richer class 



of memory effects. 

Earlier work has reported subharmonic response in the 
magnetization of ferromagnetic materials subjected to os- 
cillating magnetic fields. ^(| ^ is a specific example of 
the general phenomenon of ferroresonance. However, this 
occurs when the magnetic material is part of an electri- 
cal circuit (e.g. as the core of an inductor) that is driven 
at finite frequency. Furthermore the intrinsic response 
of the ferromagnetic material is expected to show return 
point memory. From a dynamical systems viewpoint, this 
is simply the observation that a driven nonlinear oscilla- 
tor can respond at a frequency different from the driving 
frequency, possibly depending on the initial conditions of 
the system. There has also been theoretical work on finite 
frequency hysteresis loops in ferromagnets. 0, 0, In 
contrast, we only consider H fields that vary sufficiently 
slowly to be adiabatic and our non-return point effect is 
intrinsic to the material. 

We consider the standard Edwards-Anderson spin 
glass Hamiltonian^, with Ising spins and nearest 
neighbor interactions in three dimensions, 

H = - JijSiSj -H^Si. (I) 

(i,3) 1 

The coupling strengths Ji j are chosen to be uniform ran- 
dom numbers between — I and I . The Ising spin variables 
Si take values ±1. Because we expect to apply this model 
to granular materials, we employ free boundary condi- 
tions. Other boundary conditions have been verified to 
yield similar results. We use single spin-flip dynamics at 
zero temperature, and vary the magnetic field H adia- 
batically. Thus as the magnetic field is changed, a spin 
only flips if it is energetically favorable for it to do so. 
Once a spin flips, it can render other spins unstable. If 
more than one spin is unstable, it is important to have 
an algorithm that chooses the order that spins flip in a 
deterministic manner, to avoid introducing randomness 
into the system through the dynamics. Accordingly we 
flip the spin that lowers the energy the most. This pro- 
cess is repeated until all the spins are stable, after which 
the magnetic field can be further varied. [T^ 

We start the system at large H with all the spins pos- 
itive and lower H to H m i n adiabatically. The field is 
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FIG. 1: Steady-state magnetization as a function of magnetic 
field for a spin-glass with 8 3 spins. It can be seen that two 
cycles of the magnetic field are needed for the magnetization 
(and in fact the full spin configuration of the system) to return 
to its initial value. 
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TABLE I: Probability that m cycles of the magnetic field are 
required for the magnetization to return to its initial value, 
for systems of sizes 4 3 , 8 3 and 16 3 . The number of systems 
considered for the three different sizes was x, y and z respec- 
tively. As L is increased, there is a greater probability of 
larger values of m. 



then raised to H max . Thereafter, the field is cycled re- 
peatedly between these two extremal values. After an 
initial transient period, the system reaches steady state. 
In this steady state we determine the number of cycles of 
H that are required for the system to repeat its configu- 
ration. This varies randomly, depending on the coupling 
strengths Jij. 

Clearly if H m in and H m ax are too large in magnitude, 
the magnetization will saturate at the extremal points 
of the H cycles and the behavior seen is trivial. On 
the other hand, if the range of H is too small, hysteretic 
effects are minimal and behavior of the kind above will be 
hard to observe. By trying various values of the extremal 
fields, we find that it is best to choose them to be about 
half of the saturation field. Consequently we work with 

Hfnax Urnin 1.4. 

Figure ^ shows an example of a system with linear di- 
mension L = 8, i.e. with 8 3 spins, where the steady state 
behavior is a two-cycle, i.e. two cycles of the magnetic 
field correspond to one cycle in the magnetization. Even 
though the gap between the two halves of the magnetiza- 
tion loop is much smaller than the width of the loop itself, 
with an oscillatory magnetic field, it would be straight- 
forward to detect the subharmonic response. Although 
more complicated cycles are often found, they are not 
shown in the figure for the sake of clarity. 

More generally, m cycles of the magnetic field can be 
required for the (steady state) magnetization to return 
to its initial value. The distribution of m for systems of 
different sizes is shown in Tabic [|] One observes that the 
larger the system size, the more likely it is to find long 
cycles. In fact for our largest size, with 16 3 spins, for 
about 1% of the systems we were unable to detect any 
periodicity even when checking for m up to 128. The fact 
that for the remaining systems, where periodicity was 
seen, the biggest value of m detected was 42, suggests 



the existence of truly aperiodic behavior|l7fl. 

Although Table [I] shows an increase in typical cycle 
size as the system is made larger, to some extent this 
is not caused by the dynamics becoming more complex, 
but because of the way that the cycle length is defined. 
For instance if the system were divisible into two inde- 
pendent halves, with a 6 and 7 cycle respectively, the 
full system would need 42 cycles to return to its original 
configuration. As the system is made larger, the num- 
ber of such subsystems could proliferate, driving up the 
apparent cycle length. To eliminate this kind of behav- 
ior, it is useful to examine the power spectrum of these 
systems. One measures the magnetization at H m in and 
Hmax, the extrema of the cycles of H. If the magne- 
tization goes through an m-cycle, these measurements 
will have a periodicity of 2m, while corresponding mea- 
surements of H would oscillate between H m i n and H rnax , 
exhibiting a periodicity of 2. The discrete Fourier trans- 
form of the magnetization measurements is taken and 
the power spectrum is calculated. In the case above, the 
spectrum will show peaks only at multiples of / = 1/6 
and / = 1/7, where / = 1 is the frequency of the H 
field, whereas if the system were truly indivisible, peaks 
at multiples of / = 1/42 would be observed. In addi- 
tion, the power spectrum is crucial to detecting these ef- 
fects experimentally, as we have discussed in connection 
with figure^ Although in an experimental situation, one 
would measure the power spectrum for the full time de- 
pendent magnetization M(t) instead of considering only 
its values at the extrema of H, the discrete Fourier trans- 
form computed here should be a good indicator of the 
experimental results. 

Figure shows the power spectrum for L = 16. This 
is an average of the spectra for all the different m-cycles, 
weighted according to their probabilities (given in Ta- 
ble Sharp peaks are seen at discrete frequencies, with 
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FIG. 2: Average power spectrum for systems of 8 3 spins. The 
peaks at / = 1 and / = 0, whose heights are approximately 
3 x 10 4 , have been suppressed. The power is concentrated at 
/ = 1/2, with smaller peaks at other fractions. 
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FIG. 3: Low frequency average power spectrum for systems 
of 4 3 and 16 3 spins. (To avoid confusion with peaks, no tick 
marks are shown on the x axis.) There is substantially more 
power at low frequencies for 16 3 spins, and the spectrum ex- 
tends down to 1// ~ 40, i.e. a 40-cycle that is indivisible. 



the dominant peak at / = 1/2. Figure [3] is an enlarged 
plot of the power spectrum at low frequencies for L = 4 
and L = 16. Although small, the low frequency power 
is clearly greater for L = 16. In addition, the peaks go 
down to lower frequencies, demonstrating that the dy- 
namics get more complex as L is increased in addition to 
the trivial growth in cycle length discussed in the previ- 
ous paragraph. 

In Figure |21 we have not plotted the two trivial peaks, 



at / = and 1. These are approximately 10 5 times 
the strength of the peak at / = 1/2. As one would ex- 
pect from Table, [fl the corresponding power spectrum 
for L = 4 has much fewer peaks. However, the strength 
of the peaks is larger: the peak at / = 1/2 is approx- 
imately 10~ 4 times the strength of the peaks at / = 
and 1. More generally, as L is varied, the strength of 
the trivial peaks will be proportional to the number of 
spins, i.e. L 3 . Although the limited range of L that we 
have been able to simulate does not allow us to determine 
the dependence on L of the other peaks, for large L we 
might expect any system to behave as a collection of rea- 
sonably independent regions, and the strength of the low 
frequency peaks to be proportional to l?l 2 . However, 
this argument is hazardous, since even a small coupling 
between different regions of the system might affect sub- 
tle details of the dynamics. On the basis of our numerics, 
we make the more conservative claim, that the strength 
of the low frequency peaks increases with L, while their 
relative strength (compared to the peaks at / = and 1) 
decreases with L. In a phase-locked experimental setup, 
the absolute rather than the relative strength of these 
peaks is important. 

We now turn to experimental considerations. 
Nanoparticles of spin glass material, such as CuMn, 
with of order 16 3 separate spins would have a linear 
dimension of a few nanometers, which could be readily 
manufactured. In fact it should also be possible to 
reduce the grain size somewhat, although as discussed 
in the previous paragraph the subharmonic content is 
reduced. In order to eliminate possible problems from 
interaction between the grains, it would advisable to 
coat or embed the grains in a neutral material to isolate 
them sufficiently from each other. With N g independent 
grains, the strength of the signal would be proportional 

tO ^fNg. 

As discussed earlier, the field should be varied over a 
range of order half the saturation field, which would be 
~ IT. The time dependence of the magnetic field should 
be sufficiently slow so as to evolve the spins adiabatically, 
and to prevent any significant eddy current heating. For 
the small grains that we consider, a frequency of the order 
of 10Hz should be adequate. 

The critical temperature T c of CuMn is of order IK; in 
order to see the effects that we have discussed, which are 
very small compared to the overall hysteresis loop, it is 
necessary to essentially eliminate thermal noise over the 
time scales we are interested in. With activated dynam- 
ics, it should be possible to achieve this if the temperature 
is of order lOmK. 

It should be possible to satisfy all these conditions 
in an experiment without too much difficulty. We note 
that recent work has found super s pin glass behavior in 
nanoparticle multilayer composites [18| and in frozen fer- 
rofluids 0| , where the effective individual spins are large 
and T c is correspondingly high. While this would render 
the phenomena we wish to observe accessible at higher 
temperatures, one should be careful because the satura- 
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tion field increases with T c . 

From a theoretical perspective, it is of interest to ask 
whether the results we have found here are generic for 
systems with hysteresis. As mentioned earlier in this pa- 
per, there are strong constraints on many ferromagnetic 
models that force them to have return point memory, 
precluding the possibility of multi-cycle hysteresis loops. 
Even with the spin glass model that we have considered 
in this paper, we find numerically that in one dimensional 
systems, when H is cycled repeatedly, the hysteresis loop 
always settles down to a one-cycle. However, steady state 
multi-cycles are seen in two dimensions, at rates compa- 
rable to those for three dimensional systems. 

From a dynamical systems approach, it would be nat- 
ural to expect a system with many interacting degrees of 
freedom to exhibit complex dynamical behavior, possibly 
even chaos. In order to verify whether this is indeed the 
source of the results reported here, we have carried out 
similar numerical simulations on the Sherrington Kirk- 
patrick model |2(| for spin glasses, as an example of a 
coupled non-linear system with several degrees of free- 
dom (the spins). We find that one needs a minimum of 
5 spins to find multi-cycle steady state magnetization. 
In this case, one finds only three cycles and no others 
aside from the trivial closed loop case. The fraction of 
systems showing three cycles for 5 spins is low, approx- 
imately 3 x 10~ 4 . However this fraction rises steeply 
with the number of spins. One has to be careful in ap- 



plying general results from dynamical systems, since the 
dynamics we consider here are of a special kind, being 
driven by energy minimization. For instance, with a 
time-independent magnetic field and any initial config- 
uration, the system must evolve towards a fixed point. 
Further work in this direction may explain the essential 
ingredients necessary to see the behavior reported in this 
paper. 

Our results for the SK model suggest that systems 
which are described by a Landau theory with five or more 
components to the order parameter might, in the correct 
parameter regime, show multi-cycle hysteresis behavior. 
This would have the advantage that the effect would show 
up in intensive quantitites. As discussed earlier in this 
paper, for the spin-glass system we have considered here 
the non-return point effects for the average magnetiza- 
tion per spin should vanish in the large size limit. |(| 

In conclusion, we have shown in this paper that hys- 
teresis loops have much richer behavior than commonly 
assumed. For Ising spin glass models with zero temper- 
ature dynamics, a slowly varying cyclical magnetic field 
produces a subharmonic response in the magnetization. 
Our estimates indicate that this effect should be observ- 
able in low temperature experiments. Broader implica- 
tions for other systems have been discussed. 

We thank David Belanger, Peter Young and Steve Ford 
for useful discussions. 
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